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Taking into account the electromagnetic and gravitational form factors, calculated from a 
new set of i-dependent GPDs, a new model including the soft and hard pomerons is build. 
In the framework of this model the qualitative description of all existing experimental 
data at yfs < 52.8 GeV, including the Coulomb range and large momentum transfers, is 
obtained with only 4 free parameters. Predictions for LHC energies are made. 

1 Introduction 

The dynamics of strong interactions finds its most complete representation in elastic scattering 
at small angles. Only in this region of interactions can we measure the basic properties that 
define the hadron structure: the total cross section, the slope of the diffraction peak and the 
parameter p(s,t). Their values are connected on the one hand with the large-scale structure 
of hadrons and on the other hand with the first principles which lead to the theorems on the 
behavior of the scattering amplitudes at asymptotic energies [TJ [2] • 

The definition of the structure of the high-energy elastic hadron-hadron scattering amplitude 
at small angles in the new superhigh energy range is a fundamental research problem. In this 
kinematical domain PQCD cannot be directly operative; however, there are numerous results 
derived in the framework of axiomatic field theory that can guide us both at the theoretical 
and experimental levels. 

The hard pomeron which is obviously present in deeply inelastic processes with the large 
£2 ~ 0.4 [3j will lead to a strong decrease of the energy at which the saturation or black-disk 
regime appears. It is not obvious how the total cross sections will grow with energy, especially 
in the energy region of the LHC. In the present work I investigate the impact of the hard 
pomeron on some of the features of elastic proton-proton scattering at LHC energies and small 
momentum transfer. 

The situation is complicated by the possible transition to the saturation regime, as the Black 
Disk Limit (BDL) will be reached at the LHC [H [5] . The effect of saturation will be a change in 
the i-dependence of B and p, which will begin for y/s = 2 to 6 TeV, and which may drastically 
change B{t) and p{i) at y/s = 14 TeV J5J E] ■ As we arc about to see, such a feature can be 
obtained in very different models. 

There are indeed many different models for the description of hadron elastic scattering at 
small angles [7] They lead to different predictions for the structure of the scattering amplitude 
at asymptotic energies, where the diffraction processes can display complicated features [8]. 
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This concerns especially the asymptotic unitarity bound connected with the Black Disk Limit 
(BDL) [19]. 

The total helicity amplitudes can be written as $i(s,t) = <I>' l (s,i) + $f m (s, t)e v ( s '*) , where 
$>i(s,t) comes from the strong interactions, <j>^ m (s,i) from the electromagnetic interactions 
and <p(s, t) is the interference phase factor between the electromagnetic and strong interactions 
[9] [TQl [IT] . For the hadron part the amplitude with spin- flip is neglected, as usual at high 
energy. 

In practice, many different partial waves with I — > oo must be summed and this leads to 
the impact parameter representation |13| converting the summation over I into an integration 
over the impact parameter b. In the impact-parameter representation the Born term of the 
scattering amplitude will be 

X(s,b) ~ J cPqe*-? F Bom (s,q 2 ) , (1) 

where t = —q 2 and dropping the kinematical factor 1 / yjs(s — 2m 2 ) and a factor s in front of 
the scattering amplitude. After unitarisation, the scattering amplitude becomes 

F(s,t) ~ J e^T{s,b) d 2 b. (2) 

The overlap function T(s,b) can be a matrix, corresponding to the scattering of different spin 
states. Unitarity of the S'-matrix, SS + < 1, requires that T(s,b) < 1. There can be different 
unitarization schemes which map x(sj^) to different regions of the unitarity circle [14] , In this 
work I used the standard eikonal unitarisation scheme which leads to the standard regime of 
saturation, i.e. the BDL [T5] : 

F( s ,6) = l-ex P Hx(s,fo)]. (3) 



2 Born amplitude in the impact-parameter representation 

In different models one can obtain various pictures of the profile function based on different rep- 
resentations of the hadron structure. In this model I suppose that the elastic hadron scattering 
amplitude can be divided in two pieces. One is proportional to the electromagnetic form factor. 
It plays the most important role at small momentum transfer. The other piece is proportional 
to the matter distribution in the hadron and is dominant at large momentum transfer. 

As in the EPSH model [Bj, I take into account the contributions of the soft and hard 
pomcrons. In this approach the nucleon-nucleon elastic scattering amplitude is proportional to 
the electromagnetic hadrons form-factors and can be approximated at small t by 

T(s,t) = [ ki (s/s ) ei e^ * ln[s,sa) + k 2 (s/s ) e2 e^ 4 ln ^/ s ^} G 2 em {t), (4) 

where k\ = 4.47 and k 2 = 0.005 are the coupling of the "soft" and "hard" pomerons, and 
6i = 0.00728, a[ = 0.3, and e 2 = 0.45, a' 2 = 0.10 are the intercepts and the slopes of the two 
pomeron trajectories. The normalization sq will be dropped below and s contains implicitly 
the phase factor exp(— in/2). I shall examine only high-energy nucleon-nucleon scattering with 
>/s > 52.8 GeV. So, the contributions of reggeons and odderon will be neglected. This model 
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only includes crossing-symmetric scattering amplitudes. Hence the differential cross sections of 
the proton-proton and proton-antiproton elastic scattering are equal. 
The assumption about the hadron form-factors leads to the amplitude 

T(s, t)Bom. = h\{F Born + F Born )G 2 m + h%{F Born + F Born )G 2 rav , (5) 

I suppose a non-linear trajectory for the pomeron and, as a first approximation, assume that 
the coupling is proportional to the gravitational form factor and that both soft and hard terms 
in the FBom(s, t) have a = at large t. 

3 Hadron form factors 

As was mentioned above, all the form factors are obtained from the GPDs of the nucleon. The 
electromagnetic form factors can be represented as first moments of GPDs 

F?(t)= [ dxH q {x,t); F$(t)= [ dx £ q (x,t), (6) 
Jo Jo 

following from the sum rules [TH [T7] . 

Recently, there were many different proposals for the t dependence of GPDs. We introduced 
a simple form for this t-dependence [18j . based on the original Gaussian form corresponding 
to that of the wave function of the hadron. It satisfies the conditions of non-factorization, 
introduced by Radyushkin, and the Burkhardt condition on the power of (1 — x) n in the 
exponential form of the i-dependence. With this simple form we obtained a good description of 
the proton electromagnetic Sachs form factors. Using the isotopic invariance we obtained good 
descriptions of the neutron Sachs form factors without changing any parameters. 

The Dirac elastic form factor can be written 

G 2 (t) = h fa e dl 1 + h fb e d * 1 + h fc e d3 ». (7) 

with hf a = 0.55, hfb = 0.25, hf a = 0.20, and di = 5.5, = 4.1, = 1.2. The exponential form 
of the form factor lets us calculate the hadron scattering amplitude in the impact parameter 
representation [6]. 

I shall use this model of GPDs to obtain the gravitational form factor of the nucleon in the 
impact-parameter representation. This form factor can be obtained from the second momentum 
of the GPDs. Taking instead of the electromagnetic current J M the energy- momentum tensor 
Tfiv together with a model of quark GPDs, one can obtain the gravitational form factor of 
fermions 

J dx x[H(x, A 2 , ± E(x, A 2 , £)] = A q (A 2 ) ± B g (A 2 ). (8) 
For £ = one has 

/ dx x[H(x,t)±£{x,t)}= A q (t)±B q (t). (9) 
Jo 

Calculations in the momentum-transfer representation show that the second moment of the 
GPDs, corresponding to the gravitional form-factor, can be represented in the dipolc form 

A(t) = L 2 /(l-t/L 2 ) 2 . (10) 
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Figure 1: da/dt at s/s = 52.8 GeV for pp elastic scattering, at small \t\ (left) and at large \t\ 
(right). 



with the parameter L 2 = 1.8 GeV 2 . For the scattering amplitude, this leads to 

A(s,b) = ^-K 3 (Lb), (11) 

where K^(Lb) is the MacDonald function of the third order. To match both parts of the 
scattering amplitude, the second part is multiplied by a smooth correction function which 
depends on the impact parameter 

m = (i + yfi + by^fi + b 2 ). (12) 



4 Description of the differential cross sections 

The model has only four free parameters, which are obtained from a fit to the experimental 
data: 

hi = 1.09 ± 0.004; h 2 = 1.57 ± 0.006; r\ = 1.57 ± 0.02; r 2 = 5.56 ± 0.06. 

I used all the existing experimental data at y/s > 52.8 GeV, including the whole Coulomb 
region and up to the largest momentum transfers experimentally accessible. In the fitting 
procedure, only statistical errors were taken into account. The systematic errors were used as 
an additional common normalization of the experimental data from a given experiment. As 
a result, one obtains J2xt/N — 3. where ./V = 924 is the number of experimental points. Of 
course, if one sums the systematic and statistical errors, the x 2 /N decreases, to 2. Note that 
the parameters are energy-independent. The energy dependence of the scattering amplitude is 
determined only by the intercepts of the soft and hard pomerons. 

In Fig. 1 the differential cross sections for proton-proton elastic scattering at y/s = 52.8 GeV 
are presented. At this energy there are experimental data at small (beginning at —t = 
0.0004 GeV 2 ) and large (up to —t = 10 GeV 2 ) momentum transfers. The model repro- 
duces both regions and provides a qualitative description of the dip region at -Is 1.4 GeV 2 , 
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Figure 2: da/dt for pp elastic scattering at small \t\, at \fs = 541 GeV (left) and °Js = 1800 GeV 
(right) 

for = 53 GeV 2 and for ^/s = 62.1 GeV 2 , although only the low- t abd high- t regions have 
been fitted to. 

Now let us examine the proton-antiproton differential cross sections. In this case at low 
momentum transfer the Coulomb-hadron interference term plays an important role and has 
the opposite sign. The model describes the experimental data well. In this case, the first part 
of the scattering amplitude determines the differential cross sections, and is dominated by the 
exchange of the soft pomeron. 

There was a significant difference between the experimental measurement of p, the ratio 
of the real part to the imaginary part of the scattering amplitude, between the UA4 and 
UA4/2 collaborations at y/s = 541 GeV. As shown in Table 1, the resulting values for p(0) 
appear inconsistent. A more careful analysis [251 126] shows that there is no contradiction 
between the measurements of UA4 and UA4/2. Now the present model gives for this energy 
p(y/s = 541GeV, t = 0) = 0.163, so, practically the same as in the previous phenomenological 
analysis. 



p(y/s = 541 GeV, t = 0) 


experiment 


experimental analysis 


global analysis 


This model 


UA4 


0.24 ±0.02 


0.19 ±0.03 [El]] 




UA4/2 


0.135 ± 0.015 


0.17 ±0.02 [[26]] 


0.163 



Table 1: Average values of p, derived with fixed total cross section (first two columns), and 
from a global analysis (last two columns). 

Now let us examine the data at higher energy, where the contribution of the hard pomeron 
will be more important. In Fig. 2 the description of the proton-antiproton scattering at y/s = 
541 GeV 2 and at yfs = 1800 GeV 2 is shown. In this case the Coulomb-hadron interference 
term leads to a large value of the real part of the scattering amplitude, which is determined by 
the contribution from the hard pomeron. The good description of the experimental data shows 
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Figure 3: da/dt for pp elastic scattering at large momentum transfer, at s/s = 630 GeV (left) 
and -fs = 1800 GeV (right) 

that the parameters of the hard pomeron correspond to the real physical situation. 

Figure 3 shows the description of the experimental data at larger momentum transfers for 
y/s = 630 GeV 2 and y/s = 1800 GeV 2 . It is clear that the model leads to a good description 
of these data. However, one must note that the fine structure of the dip is not reproduced by 
the model in this case. The model shows only an essential change of the slope in this region. 

Saturation of the profile function will surely control the behavior of atot at higher energies 
and will result in a significant decrease of the LHC cross section. For the last LHC energy 
sfs = 14 TeV the model predicts atot — 146 mb and p(0) = 0.235. This result comes from the 
contribution of the hard pomeron and from the strong saturation from the black disk limit. 

5 Conclusion 

In the presence of a hard Pomeron , saturation effects can change the behavior of the cross 
sections already at LHC energies. A new model, taking into account the contributions of the soft 
and hard pomerons and using form factors calculated from the GPDs, successfully describes 
all the existing experimental data on elastic proton-proton and proton-antiproton scattering 
at ^/s > 52.8 GeV, including the Coulomb- hadron interference region, the dip region, and 
the large-momentum-transfer region. The behavior of the differential cross section at small t 
is determined by the electromagnetic form factors, and at large t by the matter distribution 
(calculated in the model from the second momentum of the GPDs) as was supposed by Mittinen 
a long time ago [27] . 

The model leads to saturation of the BDL in the TeV region of energy As a result the 
parameters of the scattering amplitude B(s,t) and of p(s,t) have a complicated dependence 
on s and t and the scattering amplitude has a non-exponential behavior at small momentum 
transfer. 

The possibility of a new behavior of p(s, t) and B(s, t) at LHC energies [5] has to be taken 
into account in the procedure extracting the value of the total cross sections by the standard 
method l28l . 



6 



EDS'09 



Acknowledgements 



I would like to thank for helpful discussions J.R. Cudell and O.V. Teryacv. I gratefully ac- 
knowledges the organization committee and R. Orava for the financial support to take a part 
in the conference and and would like to thank the FRNS and University of Liege where part of 
this work was done. 

References 

[1] A.Martin, F. Cheung, Analytic properties and bounds of the 

[2] S.M.Roy, Phys.Rcp. C 5 (1972) 125 . 

[3] A. Donnachic and P. V. Landshoff, Phys. Lett. B 595 (2004) 393 [arXiv:hep-ph/040208l| . 

[4] O. V. Selyugin and J. R. Cudell, Czech. J. Phys. 54 (2004) B161 [arXiv:hep-ph/0301048| . 

[5] J. R. Cudell and O. V. Selyugin, Phys. Lett. B 662 (2008) 417 [arXiv:hep-ph/0612046| . 

[6] J. R. Cudell and O. V. Selyugin, Phys. Lett. B 662 (2008) 417 [arXiv:hep-ph/0612046| . 

[7] R. Fiore, L. Jenkovszky, R. Orava, E. Predazzi, A. Prokudin, O. Selyugin, Mod. Phys., A24: 2551-2559 
(2009) 

[8] J. R. Cudell and O. V. Selyugin, Czech. J. Phys. 54 (2004) A441 [arXiv:hep-ph/0309194| . 

[9] O.V. Selyugin, Mod. Phys. Lett. A9 (1994) 1207. 

[10] O.V. Selyugin, Mod. Phys. Lett. A14, 223 (1999). 

[11] O. V. Selyugin, Phys. Rev. D 60 (1999) 074028 

[12] G. Cohen- Tannoudji, Ph. Salin, and A. Morel, Nuov. Cim. 55 A (1968) 412 

[13] E. Predazzi, Ann. of Phys., 36 (1966) 228, 250 

[14] J. R. Cudell, E. Predazzi and O. V. Selyugin, Phys. Rev. D 79 (2009) 034033 | arXiv:0812.0735l [hep-ph]] . 

[15] O. V. Selyugin, talk at the Internatianal Bogoliubov conference "Problems of theoretical and mathematical 
physics", August 21-27, 2009, Moscow-Dubna, Russia. 

[16] X.D. Ji, Phys. Lett. 78 , (1997) 610; Phys. Rev D 55 (1997) 7114. 

[17] Radyushkin, A.V., Phys. Rev. D 56, 5524 (1997) . 

[18] O.V. Selyugin, O.V. Teryaev, Phys. Rev., D (2008). 

[19] O. V. Selyugin, J. R. Cudell and E. Predazzi, Eur. Phys. J. ST 162 (2008) 37 |arXiv:0712^62T1 [hep-phll. 
[20] UA4/2 Coll., C. Augier et al, Phys. Lett. B 316, (1993) 448. 

[21] O. V. Selyugin and J. R. Cudell, Nucl. Phys. Proc. Suppl. 146 (2005) 185 [arXiv:hep-ph/0412338] . 
[22] L. Van Hove, Rev. Mod. Phys. 36 (1964) 655. 

[23] J. R. Cudell, A. Le ngyel, E. Martynov and O. V. Selyugin, Nucl. Phys. A 755 (2005) 587 
[arXiv:hep-ph/0501288] . 

[24] J. R. Cudell, E. Ma rtynov, O. V. Selyugin and A. Lengyel, Phys. Lett. B 587 (2004) 78 
[arXiv:hep-ph/0310198] . 

[25] O. Selyugin, Sov. J. Nucl. Phys. 55, 466 (1992) 

[26] O. Selyugin, Phys. Lett. B333, 245 (1994) 

[27] Mittinen 

[28] J. R. Cudell and O. V. Selyugin, Phys. Rev. Lett. 102 (2009) 032003 |arXiv:0812l892l [hep-ph]]. 



EDS'09 



7 



